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1.  Introduction 

It  is  important  to  be  able  to  assess  the  reliability  of  a  complex  system,  based  on 
knowledge  concerning  the  reliabilities  of  its  individual  components.  Such  problems 
arise  with  increasing  frequency  in  the  analysis  of  telecommunication,  distribution  and 
transportation  systems,  all  of  which  can  be  modeled  as  networks  of  various  sites  (or 
nodes)  interconnected  by  edges.  In  a  typical  situation,  edges  of  the  network  (which 
may  be  directed  or  undirected)  are  assumed  to  fail  in  a  statistically  independent 
fashion  with  known  probabilities.  For  such  networks,  a  variety  of  probabilistic 
measures  of  system  performance  have  been  considered:  namely,  the  probability  that 
two  given  nodes  can  communicate  [1 1 ,13,33],  that  a  given  node  can  communicate  with 
a  specified  portion  of  the  network  [6,25],  or  that  every  pair  of  nodes  can  communicate 
[8,26,33], 

Numerous  algorithms  have  been  proposed  for  calculating  these  measures  of 
network  reliability.  One  class  of  methods  is  based  on  the  idea  of  a  path,  a  minimal  set 
of  edges  whose  operation  ensures  that  the  system  functions.  In  this  approach,  paths 
must  first  be  enumerated  and  then  combined  either  by  applying  the  inclusion-exclusion 
principle  or  by  effecting  a  partition  into  mutually  disjoint  events  [2,13-15,1 7,27].  An 
alternative  approach  uses  instead  the  enumeration  and  combination  of  cutsets, 
minimal  sets  of  edges  whose  failure  ensures  that  the  system  cannot  function 
[1,11,12,17,22]. 

A  different  approach,  called  pivotal  decomposition,  considers  the  two  states 


(operative  or  failed)  assumed  by  a  given  edge  in  the  network  [2,7,16,26].  Based  on 
these  two  states,  the  given  network  decomposes  into  two  smaller  subnetworks.  This 
process  is  repeatedly  applied  until  all  resulting  subnetworks  are  solvable. 
Satyanarayana  and  Chang  [26]  have  formulated  and  implemented  an  optimal  edge 
selection  strategy  that  achieves  the  minimum  possible  computational  effort  in  the  case 
of  undirected  networks. 

Despite  the  considerable  effort  spent  in  developing  algorithms  for  the  calculation 
of  network  reliability,  the  effective  computation  of  such  probabilistic  measures  for 
general  networks  remains  elusive.  Namely,  all  known  general  procedures  exhibit  a 
worst-case  behavior  that  is  exponential  in  the  size  of  the  network.  Since  no 
polynomially-bounded  algorithm  has  yet  been  devised,  the  exact  computation  of 
network  reliability  has  been  confined  to  networks  of  rather  small  size.  In  fact,  most 
reliability  problems  of  any  substance  are  now  known  to  be  NP-hard  or  #P-complete 
[5.20-22,24,32],  making  it  extremely  unlikely  that  any  polynomially-bounded  algorithm 
will  ever  be  found.  As  a  result,  recent  research  has  concentrated  on  developing 
algorithms  adapted  to  special  network  structures  (where  polynomial-time  algorithms 
are  possible),  or  has  pursued  a  simulation  approach  [10]. 

A  number  of  special  classes  of  undirected  networks  have  recently  been  analyzed 
with  success.  Polynomial-time  algorithms  are  now  available  for  calculating  certain 
reliability  measures  in  series-parallel  [28],  inner-four-cycle-free  [18],  and  cube-free  [19] 
planar  graphs.  Provan  [20]  has  shown,  however,  that  the  problem  of  determining 
source-to-terminal  reliability  remains  #P-complete  for  the  general  class  of  planar 
graphs.  In  order  to  analyze  more  complex  networks  topologies,  the  idea  of  pivotal 
decomposition  together  with  polygon-to-chain  reductions  [34]  can  be  used  to 
decompose  the  original  problem  into  smaller  subproblems. 

Similar  results  and  tools  are  not  as  readily  available  in  the  case  of  directed 
networks.  The  only  significant  types  of  directed  networks  that  are  known  to  admit  a 
polynomial-time  algorithm  are  "basically-series-parallel"  networks  [3,4].  For  general 


directed  networks,  the  algebraic  approach  developed  by  Shier  and  Whited  [30,31] 
yields  algorithms  for  solving,  either  approximately  or  exactly,  some  relatively  complex 
source-to-terminal  problems,  although  the  computation  time  grows  relentlessly  with  the 
size  of  the  network. 

This  paper  investigates  a  general  algorithm  for  calculating  source-to-terminal 
reliability  in  certain  well-structured  classes  of  networks,  exemplified  by  planar  graphs. 
The  emphasis  is  on  identifying  an  underlying  lattice  structure  that  captures  certain 
algorithmically  desirable  features  of  such  networks.  The  motivation  for  the  present 
work  derives  from  several  recent  results  concerning  cutset-based  algorithms  for 
network  reliability.  Specifically,  Buzacott  [9]  discusses  how  reliability  formulas  using 
the  disjoint  products  form  can  be  simplified  by  making  use  of  a  certain  ordering  of 
cutsets  in  the  network.  Shanthikumar  [29]  has  also  exploited  ordering  properties  of 
cutsets  to  obtain  an  efficient  recursive  algorithm  to  calculate  an  upper  bound  on 
source-to-terminal  reliability.  Most  recently.  Provan  and  Ball  [22]  have  described  a 
pseudo-polynomial  algorithm  for  computing  the  reliability  of  an  arbitrary 
source-to-terminal  network  (either  directed  or  undirected).  Namely,  an  algorithm  is 
provided  with  worst-case  complexity  that  is  polynomial  in  the  number  of  cutsets 
separating  the  source  s  and  terminal  t.  Provan  and  Ball  further  show  that  for  directed 
networks  no  algorithm  polynomial  in  the  number  of  (s,t)  paths  can  exist  unless  P  =  NP. 

For  networks  of  reasonable  size  and  sparsity,  the  collection  of  (s.t)  paths  can  be 
quickly  enumerated  using  a  modified  depth-first  search  of  the  network  [23].  The 
generation  of  all  (s.t)  cutsets  is  not  as  straightforward,  particularly  in  the  case  of 
directed  networks.  Thus,  while  a  pseudo-polynomial  algorithm  exists  relative  to  the 
number  of  cutsets,  the  effort  required  to  generate  such  cutsets  may  be  strikingly  larger 
than  that  required  to  generate  the  corresponding  paths.  As  an  example,  consider  the 
dodecahedron  graph  given  in  Figure  1 ,  having  20  nodes  and  30  edges.  This  example 
has  7376  (s,t)  cutsets  but  only  780  (s.t)  paths.  Thus,  in  this  instance  not  only  are  the 
cutsets  more  difficult  to  generate  than  the  paths,  but  there  are  substantially  more 
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cutsets  than  paths  involved. 

The  emphasis  here  is  to  crystallize  those  features  of  cutset-based  algorithms  that 
make  them  theoretically  attractive  (pseudo-polynomial  behavior),  while  at  the  same 
time  broadening  the  scope  from  cutsets  to  other,  more  easily  computed,  entities.  We 
present  a  general  framework  and  a  general  recursive  algorithm  defined  on  certain 
collections  of  objects  (e.g.,  cutsets  and  paths)  that  admit  a  suitable  partial  ordering. 
This  algorithm  will  run  in  time  polynomial  in  the  number  of  such  objects. 
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2.  Basic  Definitions  and  Concepts 

Let  E  =  be  a  set  of  components,  subject  to  failure,  and  let  X  =  {X, . Xr}  be 

a  collection  of  subsets  of  E.  Components  in  E  have  two  states:  active  and  inactive.  A 
set  Xj  c  E  is  called  active  if  all  its  components  are  active.  We  suppose  that  (X,>)  is  a 
partial  ordering  having  the  lattice  property:  namely,  any  two  Xjt  Xj  e  X  have  a  unique 
least  upper  bound  Xj  v  Xj  and  a  unique  greatest  lower  bound  X(  a  Xj.  There  are  two 
additional  requirements  imposed  here. 

Closure-  If  Xj  is  active  and  Xj  is  active,  then  Xf  v  Xj  and  Xj  a  Xj  are  active.  (1 ) 

Connectivity-  If  e  e  Xj  and  e  e  Xj  then  e  e  [Xj,  Xj]  =  { Z  e  X:  Xj  <  Z  <  Xj}.  (2) 

The  event  {Xj  is  active}  will  be  denoted  by  Aj.  We  will  be  interested  in  calculating  the 
quantity 

Q(X)  =  Pr(At  u  A2  u ...  v  Ar). 


Example  1.  Consider  the  source-to-terminal  undirected  network  G  shown  in  Figure  2, 
with  components  (edges)  E  =  {e^, 63,64,65}.  Here  the  Xj's  will  be  the  (s,t)  cutsets  in 
the  network  defined  by  the  following  edge/node  sets: 


X1={e1,e2}  V,  =  {s} 

X2  =  {©1  >©3.65}  V2  =  {s.b} 

X3  =  {62.63,64}  V3  =  {s,a} 

X4  =  {64,85}  V4  =  {s,a,b} 


That  is,  the  cutset  X,  consists  of  all  edges  in  G  =  (V,E)  joining  nodes  of  V(  with  nodes  of 
V-Vj.  The  partial  order  >  is  defined  by 
Xj  >  Xj  <=>  V;  2  Vj  . 

The  associated  Hasse  diagram  for  this  partial  order  is  shown  in  Figure  3.  It  is  easy  to 


verify  that  requirements  (1)  and  (2)  both  hold.  By  considering  an  "active"  component  to 
mean  a  "failed"  component,  the  event  Aj  signifies  that  all  components  in  Xj  fail  and  n(X) 
is  just  the  network  unreliability  1-Rsl(G):  i.e.,  the  probability  that  a  message  sent  from  s 
will  be  unable  to  reach  t. 

Example  2.  Again  consider  the  network  in  Figure  2,  with  the  Xj's  now  denoting  the  (s.t) 
paths: 

Xi  =  {02'%} 

X2  =  {61.63.65} 

X3  =  {62.63,64} 

X4  =  {6^64} 

We  define  Xj  >  Xj  to  mean  that  path  Xf  is  geometrically  "above"  path  Xj,  yielding  the 
same  Hasse  diagram  shown  in  Figure  3.  Again,  requirements  (1)  and  (2)  are  seen  to 
hold.  Here,  the  association  of  "active"  with  "functioning"  allows  Aj  to  be  interpreted  as 
the  event  that  all  components  in  Xs  are  functioning,  whereupon  Q(X)  is  the 
source-to-terminal  reliability  Rst(G). 

As  the  above  examples  illustrate,  we  can  define  a  partial  ordering  on  the  (s.t) 
cutsets  in  an  arbitrary  network  or  on  the  (s.t)  paths  in  an  (s.t)-planar  network.  In  both 
cases,  stipulations  (1)  and  (2)  are  guaranteed  hold.  Other  (non-network)  examples  of 
coherent  systems  [7]  also  exist  which  admit  an  ordering  such  that  (1)  and  (2)  are 
satisfied  with  respect  to  the  system's  "minpaths"  or  "mincuts".  The  following  example 
illustrates  one  such  case. 

Example  3.  Consider  a  system  with  components  {e^,. ...eg}.  The  system  functions  if 
and  only  if  all  components  function  in  at  least  one  of  the  following  subsets  ("minpaths"): 
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X,  =  {©.j  i©2>®3} 

X2  —  {02’®3'®4} 

X3  =  {©3.64.65} 

X4  =  {64.65.65} 

If  the  (total)  order  >  is  defined  by  X4  s  X3  £  X2  £  X,  then  requirements  (1 )  and  (2)  will 
hold.  In  this  case,  the  Hasse  diagram  for  the  partial  order  is  simply  a  "chain,"  and  the 
sets  Xj  do  not  correspond  to  the  paths  or  cutsets  of  any  network. 


§ 
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3.  The  Recursive  Algorithm 

In  what  follows,  it  will  be  assumed  without  loss  of  generality  that  the  elements  of 
the  partial  order  have  been  topologically  ordered :  that  is, 

Xj  <  Xj  =>  i  <  j  .  (3) 

We  are  interested  in  calculating  Q(X)  =  Pr(A),  where  A  =  A1uA2u..,u  Ar  The  events 
Aj  are  not  disjoint,  however.  Accordingly,  we  define  the  events  Fj  (which  will  be 
disjoint)  by 

Fj  =  {Xj  is  the  "lowest"  active  set  in  X). 

These  events  are  well  defined  in  view  of  (1).  Namely,  if  Xj^  Xj . Xik  are  active  sets 

then  so  is  X:  a  X:  a  ...  a  X=  ,  and  thus  X=  aX.  a...  aX=  =  X:  for  some  j.  Notice  that 
X:  aX:  a...  aX  <X=  for  all  m,  and  so  X:  ^  X:  for  all  active  sets  X=  .  Thus  X,  is 

'1  '2  'k  'm  'j  'm  'm  'j 

indeed  the  lowest  active  set  in  X. 

Our  general  development  will  parallel  that  of  Provan  and  Ball  [22],  which  was 
originally  stated  in  the  specific  context  of  (s,t)  cutsets.  First,  it  is  easy  to  establish  the 
following  two  properties  of  the  Fj's. 

Property  1 .  FjnFj^Ofori  *j. 

Proof:  If  Fj  and  Fj  both  occur  for  i  ^  j  then  Xj  and  Xj  are  active,  so  Xj  a  Xj  is  active  by  (1 ). 
Since  Xj  a  Xj  <  X;  and  Xj  a  Xj  <  Xj  then  Xj  =  Xj  a  Xj  =  Xj,  a  contradiction.  ♦ 

Property  2.  A  =  F,  u  F2  u ...  u  Fr . 

Proof:  Suppose  A  occurs  with  X^  being  active.  Then  X(l  a  ...  a  Xik  =  X,.  is  active 
by  (1 )  and  so  Fj.  occurs.  Conversely,  if  Fj  occurs  then  Xj  is  active  and  so  A  occurs.  ♦ 

As  a  result  of  the  above  properties 
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Q(X)  =  Pr(A)  =  X  Pr(R) , 


(4) 


i-i 


and  the  problem  reduces  to  that  of  calculating  the  Pr(Fj).  A  general  recursion  involving 
these  quantities  can  be  easily  derived  as  follows. 

Aj  =  U  (FjO  Aj) 


i-i 


=  [X|VXj(^Ai)]-,i 

=  [  LJ  Fjn{X:-X.  is  active}]  u  F 
LX1<Xi  J  1 


(5) 


The  first  equality  follows  from  Property  2.  The  second  follows  since  if  Fj  n  Aj 


occurs  then  Xj  a  Xj  =  Xj,  whence  Xj  <  Xj.  Because  all  unions  in  (5)  involve  disjoint 
events, 


Pr(Aj)  =[  £  Pr(Fj )  Pr(X.  -  X.  is  active  |  F.)]  +  Pr(F) . 
X:  ^  X: 


Now  by  requirement  (2),  e  €  Xj  -  X(  cannot  be  an  element  of  any  Xk  <  Xj.  This  means 
that  the  event  Fjt  which  requires  Xj  to  be  active  but  all  Xk  <  Xj  to  be  inactive,  is 
independent  of  the  event  {Xj  -  Xj  is  active},  whence 


Pr(Aj)  =  [  £  Pr(Fj)ajj]  +  P(Fj)  , 


Xj<Xj 


(6) 


where 


ctjj  =  Pr(Xj  -  Xj  is  active) 


=  I!  Pr(e  is  active) , 


ee  Xj-X; 


assuming  the  independence  of  component  failures.  Rearranging  (6)  gives  the 
recursion  (see  [22]) 


Pr(F )  =  Pr(A )  -  X  Pr(F|)  <% 


(7) 


X^X, 


The  quantities  A,  and  ay  are  readily  computed  and  thus  Pr^),  Pr(F2) . Pr(Fr)  can  be 

found  in  turn  using  (7).  The  quantity  £2(X)  can  then  be  determined  via  (4).  The 
worst-case  complexity  of  the  entire  procedure  is  0(mx)  <  0(|X|2),  where  mx  indicates 
the  number  of  arcs  in  the  full  partially  ordered  set  (X,>).  Thus  the  above  procedure  is 
pseudo-polynomial:  namely,  its  running  time  is  bounded  by  a  polynomial  in  the 
number  of  elements  in  the  partial  order. 

Example  4.  Consider  the  source-to-terminal  network  given  in  the  upper  portion  of 
Figure  4.  The  ordering  of  paths  from  "top”  to  "bottom"  gives  the  Hasse  diagram  shown 
in  the  lower  portion  of  the  figure.  By  using  the  shorthand  notation  =  PrfAj),  f,  =  Pr(Fs) 
and  denoting  the  edge  reliabilities  as  pk,  relation  (7)  yields: 

f i  =  ai  ~  P2P7  ■ 

f2  =  a2  -  f,  (p5p6)  =  P2P5P6  ■  P2P5P6P7  ■ 

f3  =  a3  -  f1(p1P3)  =  P1P3P7  *  P1P2P3P7  • 

7 

and  so  forth,  yielding  Rst(G)=  ^(X)=  X  V 


4.  An  Alternative  Formulation 


Equation  (6),  rewritten  with  our  shorthand  notation,  becomes 

—  X  f  j  (Xjj  t  (Xjj  =  1  . 

1  Xj  <  Xj  1  u 

We  can  now  use  Moebius  inversion  over  the  partially  ordered  set  (X  >)  to  obtain 

<l  >  I  a,4  (8 

X,<  Xj 

where  Py  is  the  alternating  sum  of  chain  values,  taken  over  all  chains  Q  in  (X  >)  joining 
X|  and  Xj.  That  is, 

{1  ,  i  =  j 

E(-l)|Q,val(Q)  ,  i*j 

Q 

where  val(Q)  =  FI  ae  and  |Q|  denotes  the  number  of  arcs  in  chain  Q.  Algebraic 

eeQ 

simplification  of  (8)  yields 

-  I<-1>IQ''°<Qj)  .  (9 

°i 

where  the  chains  Qj  extend  downward  from  Xj  in  (X>)  and 
o(Qj)  =  ri{Pr(e):  e  e  Xj  for  some  Xf  e  Qj} 

gives  the  joint  probability  of  all  edges  appearing  in  some  Xf  along  the  chain. 
Analogously,  when  Q  ranges  over  all  chains  in  (X  >)  we  obtain 
n-  I Q 1 

Q(X)=I(-1)  a(Q)  .  (10 

Q 

As  an  illustration,  consider  the  bridge  network  in  Figure  2  with  the  associated 
partially  ordered  set  (based  on  paths)  given  in  Figure  5.  Then,  using  equation  (9) 

yields,  for  example 


U  ~  P1P4  *  (P 1  P3p4p5"*‘P  1  P2p4P5+P  1 P2P3P4)  +  ^PlP2P3P4P5> 


where  first  term  corresponds  to  the  chain  of  length  0,  the  next  (negated)  terms 
correspond  to  the  three  chains  of  length  1 ,  and  the  final  two  terms  correspond  to  the 
chains  of  length  2.  The  final  expansion  for  fl(X)  resulting  from  (10)  has  1 1  terms, 
corresponding  to  the  1 1  chains  in  (X  >).  This  number  is  to  be  compared  with  the 
24  -  1  =  15  potentially  occurring  terms  appearing  in  the  inclusion/exclusion  formula  [2] 
Thus  equation  (10)  captures  some  of  the  cancellation  found  in  the  topological 
expansion  formula  of  Satyanarayana  and  Prabhakar  [27]. 


5.  Computational  Results 


Promising  computational  results  have  been  obtained  using  the  recursive 
algorithm  described  in  Section  3.  We  present  in  this  section  a  fairly  complete  analysis 
of  two  particular  networks  (of  moderate  complexity)  that  have  thus  far  resisted  any 
thorough  analysis. 

The  first  network,  with  14  nodes  and  25  edges,  is  shown  in  Figure  6.  This 
network  (suggested  by  J.  S.  Provan)  is  known  to  be  "self-dual, K  having  426  paths 
joining  the  specified  s  and  t  as  well  as  426  cutsets  separating  s  and  t.  Since  the 
network  is  (s.t)-planar,  the  paths  can  be  ordered  as  described  in  Section  2  and  the 
recursive  algorithm  can  then  be  applied.  In  order  to  present  the  results  succinctly, 
Table  1(a)  presents  the  coefficients  of  the  reliability  polynomial  when  expressed  as  a 
function  Rst(p)  of  the  common  edge  reliability  pk  =  p.  Table  1(b)  lists  values  of  this 
polynomial  at  selected  values  of  p  throughout  the  range  [0,1].  As  expected,  Rst(0)  =  0, 
Rst(1)  =  1 ,  and  Rst(0.5)  =  0.5,  with  the  latter  value  occurring  because  the  network  is 
self-dual.  Computation  of  the  relevant  paths  and  the  polynomial  Rst(p)  required 
approximately  2.7  seconds  using  the  IBM  3081 -K  computer  at  Clemson  University. 

A  more  complex  network  is  the  dodecahedron,  shown  in  Figure  1 .  Since  this 
network  is  not  (s.t)-planar,  the  recursive  formula  was  applied  relative  to  the  ordering 
produced  using  the  (s,t)  cutsets.  Knowledge  of  the  planarity  of  the  network  was, 
however,  exploited  in  generating  and  processing  the  7376  cutsets.  Table  2(a) 
presents  the  coefficients  of  the  unreliability  polynomial  U8t  =  1  -  Rst,  now  expressed  for 
convenience  as  a  function  of  the  common  edge  failure  probability  qk  =  1  -  pk  =  q. 
Evaluation  of  U8l(q)  at  selected  values  of  q  is  given  in  Table  2(b).  The  total  time 
required  to  generate  the  cutsets  and  produce  Usl(q)  was  under  four  and  a  half  minutes 
of  CPU  time.  It  should  be  noted  that  the  only  previous  calculations  available  for  this 
example  required  several  hours  of  CPU  time  (on  a  similar  mainframe)  in  order  to 
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J  obtain  the  value  of  U8t(q)  at  two  selected  points  (q  =  0.1 ,  q  =  0.5). 
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Table  1.  Reliability  Polynomial  Rst(p)  for  Network  in  Figure  6 


(a)  Coefficients  of  Rst(p): 


Rst(P)  = 


4  p4  +  18  p5  +  30  p6  -  38  p7  -  188  p8  - 154  p9  +  128  p10  +  1456  p11 
+  658  p12  -  3252  p13  -  4952  p14  - 1 21 2  p15  +  37621  p16  -  1 082  p17 
-  213712  p1 8  +  483698  p19  -  565880  p20  +  414334  p21  -  198122  p22 
+  60628  p23  -  1 0850  p24  +  868  p25 


(b)  Evaluation  of  Rst(p)  at  selected  points: 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 


Rst(p)  :  .0006  .0131  .0774  .2428  .5000  .7573  .9226  .9869  .9994 
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Table  2.  Unreliability  Polynomial  Ust(q)  for  Network  in  Figure  1 


(a)  Coefficients  of  Usl(q): 


Ust(q)  = 


2  q3  +  6  q4  +  18  q5  +  89  q6  +  21 6  q7  -  561  q8  -  5522  q9  -  4935  q10 
+  56046  q"  +  133630  q12- 624126  q13  -  1223316  q14 
+  8168604  q15  -  4322481  q16  -  63482286  q17  +  256043573  q18 

*  560270502  q19  +  848614479  q20  -  961097704q21  +  841480923  q22 

*  577634706  q23  +  31 1 587945  q24  -  1 31 0921 76  q25  +  42244830  q26 

*  10090328  q27  +  1685643  q28  -  1 76016  q29  +  8656  q30 


(b)  Evaluation  of  U^q)  at  selected  points: 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 


U„(q)  :  .0029  .0362  .1685  .4272  .  7097  .8975  .9767  .9973  .9999 
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